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Computer tomography method for a periodically moving object 



PCT/IB2004/001061 



The invention relates to a computer tomography method in which a 
periodically moving object, in particular an organ of the body, which is located in an 
examination area is irradiated by a cone-shaped beam cluster. The invention also relates to a 
computer tomography scanner for carrying out the method and to a computer program for 
5 controlling the computer tomography scanner. 

Within the context of this invention, the expression "periodic movement" 
refers to movements in which a series of object states are repeatedly assumed by the object, 
always in the same sequence, during a measurement. An object state may in this case be 
defined, for example, by the position of the object in the examination area and the shape of 

10 the object or by a specific measured signal which is allocated to the object state. This 
measxured signal may, for example, be the signal of an electrocardiograph if the moving 
object is a heart. If an electrocardiograph supplies the same measured signals at different 
points in time, then it is assumed that the object states at these points in time are the same. 
Within the context of the invention, however, the term "periodic movements" also means 

15 movements which could be referred to as "quasi-periodic" and in which object states that are 
repeatedly assumed by the object are not exactly identical but are substantially identical. The 
same appUes to the time duration between two substantially identical, repeating object states, 
which time duration need only be substantially constant during the measurement in order to 
be able to refer to the corresponding movement as "periodic" within the context of this 

20 invention. Different period time durations are "substantially identical" if their time difference 
is small compared to the periodic time durations for the respective application. Accordingly, 
two object states are substantially identical if their difference is small compared to the 
differences which the object passes through within an entire period. Therefore, for example, 
movements of pulsating technical or biological objects such as organs of the body or veins 

25 may be referred to as "quasi-periodic" and hence also as "periodic". Moreover, the term 
"periodic movement" also encompasses attenuated oscillations or movements in which the 
object repeatedly assumes only some, that is to say not all, states at substantially identical 
time intervals. 
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In known methods of the type mentioned above, the spatial profile of the 
absorption or of the attenuation of the radiation in the periodically moving object is 
reconstructed jfrom measured values acquired using a detector unit. In this case, the 
movement of the object leads to the measured values containing information firom different 
object states, and this leads to movement artefacts in the reconstructed data record. 

Therefore, in known reconstruction methods use is only made of measured 
values which have been acquired in time ranges in which the object has moved relatively 
little, that is to say in which the object states change only slightly. In one of these methods, 
the beams firstly undergo parallel rebinning in order to reduce the computational complexity. 
In parallel rebinning, the cone-shaped beam cluster for each radiation source position is split 
on the trajectory into beam fans, the beams of which in each case lie in a plane which runs 
parallel to the central axis of the circular or helix-like trajectory. These beam fans are then 
divided into groups such that each group contains only beam fans that are parallel to one 
another. The beam fans of a group are referred to as a projection. A point in time is allocated 
to each projection, said point in time being the same as the point in time at which the 
measured values corresponding to the central beam fan of the projection were acquired. It is 
assumed that all measured values of a projection were acquired at this point in time. 

In the reconstruction of the absorption distribution firom the rebinned beams or 
measured values, use is only made of projections which were acquired during predefinable 
time ranges in which the movement of the object is relatively small. Since - apart fi:om in 
respect of the central beam fan - the point in time of a projection does not coincide with the 
actual acquisition time of a measured value of this projection, for the reconstruction use is 
firequently made of measured values whose beams do indeed come into contact with the 
object but were acquired at points in time which lie outside the predefined time ranges. This 
leads to movement artefacts which increase as the distance of the object firom the central 
beam fan increases. 

It is an object of the present invention to specify a method in which these 
movement artefacts are less pronounced. This object is achieved according to the invention 
by a computer tomography method comprising the steps 

using a radiation source (S) to generate a cone-shaped beam cluster (4) which 
passes through an examination area (13) and a periodically moving object which is located in 
the examination area (13), 

producing a relative movement between the radiation source (S) on the one 
hand and the object located in the examination area (1 3) on the other hand, where a 
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trajectory, along which the radiation source moves relative to the object, runs on an 
imaginary cylindrical surface that surrounds the object, 

using a detector unit (16) to acquire measured values which depend on the 
intensity in the beam cluster (4) on the other side of the object, during the relative movement, 
5 - recording the periodic movement of the object during the acquisition, 

reconstructing a spatial distribution of the absorption of tihie periodically 
moving object j&om the measured values with the aid of the recorded periodic movement of 
the object, comprising the steps: 

a) determining the spatial area taken up by the object in the examination area 
10 (13), 

b) subjecting the measured values to parallel rebinning in order to form a nmnber 
of groups, where the beams corresponding to the measured values of each group form beam 
fans (41 . . .45) which lie in planes that are parallel to one another and to the axis of rotation, 

c) determining for each group a measured value whose beam irradiates the spatial 
1 5 area taken up by the object, and allocating to the respective group the point in time at which 

this measured value was acquked, 

d) determining those groups whose points in time, allocated in step c), lie within 
periodic, predefined time ranges (Hi), 

e) reconstructing the absorption distribution in the object from the measured 
20 values belonging to the groups determined in step d). 

In contrast to known methods, in the invention there is determined for each 
group or projection a measured value whose beam irradiates the object, and the point in time 
at which this measured value was acquired is allocated to the projection. This leads to the 
time intervals between the point in time allocated to a projection and the acquisition times of 
25 said measured values of this projection, which irradiate the object, being smaller than in 
known methods. This reduces the incorrect use of measured values which lie outside the 
predefined time ranges. Since the time ranges are generally selected such that they correlate 
with phases of the object movement which do not contain much movement, movement 
artefacts are reduced. 

30 Claim 2 describes one refinement in which the spatial area taken up by the 

object in the examination area is determined with a low computational complexity by 
reconstmction and segmentation of the object. 

In Claim 3, for each group there is determined the measured value whose 
beam runs through the geometric center of the object. The point in time at which this 
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measured value was acquired is allocated to the respective group. This leads to a further 

reduction in the incorrect use of measured values which lie outside the predefinable time 

range, and this further reduces the movement artefacts. 

Claim 4 shows one refinement in which the periodically moving object is a 
5 heart and in which reconstructed images of a high quality can be produced inter alia with the 

aid of an electrocardiograph. 

In claim 5, the predefinable time ranges are selected such that the object 

moves less in these time ranges tiian in other time ranges, and this leads to a further reduction 

of the movement artefacts in the reconstructed images. 
10 Claim 6 defines one refinement in which the reconstruction is carried out 

using a filtered back-projection, and this leads to a good image quality of the reconstructed 

object while entailing a low computational complexity. 

In claim 7 the radiation source moves relative to the examination area on a 

helix-like or circular trajectory, and this leads to reconstruction results with a good image 
15 quality. 

A computer tomography scanner for carrying out the method is described in 

claim 8. 

Claim 9 defines a computer program for controlling a computer tomography 
scanner as claimed in claim 8. 



The invention will be further described with reference to examples of 
embodiments shown in the drawings to which, however, the invention is not restricted. 

Fig. 1 shows a computer tomography scanner which can be used to carry out 
25 the method according to the invention. 

Fig. 2 shows a flowchart of the method according to the invention. 

Fig. 3 shows a beam cluster of parallel rebinned beam fans. 

Fig. 4 shows a flowchart for a back-projection. 

Fig. 5 shows an electrocardiogram- 
30 Fig. 6 shows a flowchart for a further back-projection. 

Fig. 7 shows the dependence of an a factor on time. 

Fig. 8 shows the dependence of a p factor on time. 
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The computer tomography seamier shown in figmre 1 comprises a gantry 1 
which can rotate about an axis of rotation 14 that runs parallel to the z direction of the 
coordinate system shown in figure 1 . For this purpose, the gantry 1 is driven by a motor 2 at a 
preferably constant but adjustable angular velocity. A radiation source S, for example an 
X-ray generator, is attached to the gantry 1. Said radiation source is provided with a 
collimator arrangement 3 which masks, jfrom the radiation generated by the radiation source 
S, a cone-shaped beam cluster 4, that is to say a beam cluster which has, both in the z 
direction and in a direction perpendicular thereto (i.e. in a plane perpendicular to the axis of 
rotation), a finite dimension other than zero. 

The beam cluster 4 penetrates a cylindrical examination area 13 in which there 
is located a periodically moving object (not shown). In this example of embodiment, said 
object is a beating heart which carries out intrinsic movements and in some circximstances 
also moves back and forth on account of the patient's breathing movements. In other 
embodiments, other periodically moving organs of the body, such as the liver or the brain, 
periodically moving parts of organs of the body or periodically moving technical objects 
could also be irradiated. 

After passing through the examination area 13, the beam clxister 4 comes into 
contact with a detector unit 16 which is attached to the gantry 1 and has a detector surface 
comprising a large number of detector elements which in this embodiment are arranged in 
rows and columns in the form of a matrix. The detector colunms preferably run parallel to the 
axis of rotation 14. The detector rows are located in planes perpendicular to the axis of 
rotation, in this embodiment on an arc of a circle aroimd the radiation source S (focus- 
centered detector). However, in other embodiments they may also be formed differently, for 
example describe an arc of a circle arovmd the axis of rotation 14 or be rectilinear. Each of 
the detector elements with which the beam cluster 4 comes into contact supplies, in each 
position of the radiation source, a measured value for a beam from the beam cluster 4. 

The spread angle of the beam cluster 4, designated amax, determines the 
diameter of the object cylinder within which the object to be examined is located during 
acquisition of the measured values. The spread angle is defmed as the angle enclosed by a 
beam lying in a plane perpendicular to the axis of rotation 14 at the edge of the beam cluster 
4 with a plane defined by the radiation source S and tiie axis of rotation 14. The examination 
area 13 or the object or the patient table may be moved parallel to the axis of rotation 14 or to 
the z axis by means of a motor 5, However, the gantry may also be moved in this direction in 
an equivalent manner. When a technical object is concerned rather than a patient, the object 
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can be rotated during an examination while the radiation source S and the detector unit 16 
remain stationary. 

With the aid of the motors 2 and 5, the radiation source S and the detector unit 
16 can describe a trajectory relative to the examination area 13, said trajectory rumiing on an 
imaginary cylindrical surface. This trajectory may, for example, run in a helix-like manner 
when both motors are operating. If, on the other hand, tlie motor 5 for advance in the 
direction of the axis of rotation 14 is idle and the motor 2 is making the gantry rotate, the 
result is a circular trajectory for the radiation source S and the detector unit 16 relative to the 
examination area 13. In this example of embodiment, only the helix-like trajectory is 
considered. 

During the acquisition of the measured values, the movement of the heart is 
recorded in a known manner by means of an electrocardiograph 8. For tihis purpose, the chest 
area of a patient is connected to the electrocardiograph 8 by way of electrodes (not shown). 
In other embodiments, in particular in the case of other moving objects, the movement of the 
object may be traced in other ways. Thus, for example, the movement information could be 
obtained from the values measured by the detector unit 16 themselves, so that there is no 
need for the movement to be recorded using an additional device, such as an 
electrocardiograph. For this, firstly a kymogram is produced from the measured values, from 
which kymogram the movement can be derived in a known manner. A detailed description of 
this method can be found in "Kymogram detection and kymogram-correlated image 
reconstruction from subsecond spiral computed tomography scans of the heart," M. 
KachekieB, D.A. Sennst, W. Maxhnoser, W.A. Kalender, Medical Physics 29(7): 1489-1503, 
2002, to which reference is hereby made. 

In this example of embodiment, it is assumed that the patient is not breathing 
during the measurement. The breathing movement can thus be disregarded. In other 
embodiments, the breathing movement could be measured, for example, using a deformable 
chest strap which is connected to a breathing movement measurement device. 

The measured values acquired by the detector imit 16 are fed to a 
reconstruction and image processing computer 10 which is connected to the detector unit 16, 
for example via contactless data communication (not shown). Moreover, flie 
electrocardiogram is transferred from the electrocardiograph 8 to the reconstruction and 
image processing computer 10. The reconstruction and image processing computer 10 
reconstructs the absorption distribution in the examination area 13 and displays it on a 
monitor 1 1 for example. The two motors 2 and 5, the reconstruction and image processing 
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computer 10, the radiation soxirce the electrocardiograph 8 and the transfer of the measiared 
values from the detector xmit 16 to the reconstruction and image processing computer 10 are 
controlled by the control imit 7. The control unit 7 furthermore controls the transfer of the 
electrocardiogram from the electrocardiograph 8 to the reconstruction and image processing 
computer 10. 

In other embodiments, the acquired measured values and the measured 
electrocardiograms can be fed for reconstruction firstly to one or more reconstruction 
computers which forward the reconstructed data e.g. via a fiberoptic cable to the image 
processing computer. 

Figure 2 shows the sequence of a measxirement and reconstruction method 
which can be carried out using the computer tomography scanner shown in figure 1. 

After initialization in step 101, the gantry rotates at an angular velocity which 
in this example of embodiment is constant. However, it may also vary, for example as a 
ftinction of time or of the radiation source position. 

In step 103, the examination area or the object or the patient table is moved 
parallel to the axis of rotation and the radiation of the radiation source S is switched on so 
that the detector xmit 16 can detect the radiation from a large number of angular positions. At 
the same time as or before the switching on of the radiation source S, the electrocardiograph 
8 is activated so that an electrocardiogram is measured at the same time. 

Thereafter, the spatial area in which the heart is located within the examination 
area is determined. For this, the absorption distribution in the examination area is firstly 
reconstructed from the measured values at a low resolution vsdthout taking the 
electrocardiogram into account. There is a low resolution for example when a volume of 20 x 
20 X 20 cm^ is represented by 64^ voxels. 

In step 105, the measured values are parallel rebinned for the pmposes of 
reconstruction. By means of the parallel rebinning, the measured values are resorted and 
reinterpolated as though they had been measured using a different radiation source (an 
expanded radiation source which is arranged on part of a helix and can emit beam fans that 
are in each case parallel to one another) and using a different detector (a flat, rectangular 
"virtual detector" including the axis of rotation 14). 

This is explained in more detail with reference to figure 3. Here, 17 designates 
the helbc-like trajectory from which the radiation source irradiates the examination area. A 
fan-shaped beam cluster 43, the beams of which nm in a plane that includes the axis of 
rotation 14, is emitted from the radiation source position Sq . It is possible to consider that the 
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cone-shaped beam cluster emitted from the radiation source at position Sq is composed of a 
large number of planar beam fans which are located in planes parallel to the axis of rotation 
14 and intersect at the radiation source position Sq . figure 3 shows just one of these beam 
fans, namely the beam fan 43. 
5 Moreover, further beam fans 41, 42, and 44, 45 are shown in figure 3, said 

beam fans being parallel to the beam fan 43 and lying in planes that are parallel to one 
another and to the axis of rotation 14. The associated radiation source positions 5 '^-i 
iS| , iS'2 are taken up by the radiation source S respectively before and after it has reached the 
radiation source position Sq . 

10 The beam fans 41 to 45 form a group and define a beam cluster 70 with a tent- 

like shape. A group of beam fans is csdled a projection. For each projection there is then 
defined a rectangular, virtual detector 160 which lies in a plane that includes the axis of 
rotation 14 and is oriented perpendicular to the parallel beam fans of a projection. The comer 
points of the virtual detector 160 are the penetration points of the beams which from the outer 

15 radiation source positions come into contact with the opposite helix section, through this 
plane. For the beam cluster 70 in figure 3, "^2 ^® outer radiation source 

positions. Detector elements that are arranged in a Cartesian manner are defined on the 
rectangular detector 160, that is to say rows and columns on which the measured values are 
reinterpolated. 

20 Subsequently, in step 107, the measured values allocated to the individual 

beams are multiplied by a weighting factor which corresponds to the cosine of the cone angle 
of the respective beam. The cone angle of a beam is the angle encompassed by this beam 
with a plane oriented perpendicular to the axis of rotation 14. If said angle is small, then the 
cosine of the angle is essentially equal to 1, so that step 107 may be omitted. 

25 In step 1 09, a one-dimensional filtering with a transmission factor that 

increases in a ramp-like manner with the spatial frequency is applied to the measured values. 
For this, use is made of respectively successive values in the direction perpendicular to the 
axis of rotation 14, that is to say along a row of the virtual detector 160. This filtering is 
carried out along each row of the virtual detector for all groups of beam fans. 

30 In other embodiments the parallel rebinning could be omitted. It is then known 

to modify the filtering since the detector unit is curved for example in the maimer of an arc 
around the radiation source or around the axis of rotation. 
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In step 1 1 1 , the filtered measiared values are then used for reconstruction of 
the absorption distribution in the exanaination area by means of a back-projection. The 
individual steps of the back-projection are shown in figure 4. 

In step 201 a voxel F(x) is defined within a predefinable field of view (FOV). 
Since the reconstruction is to take place with a low resolution, the amount of voxels may be 
for example 64^ and the FOV may be 20 x 20 x 20 cm^. In step 203, a projection, that is to 
say a group of beam fans, is then selected which has not yet been used to reconstruct the 
voxel K(x) . If no beam of the projection runs centrally tibrough the voxel K(x) , then the 
point at which a central beam would have come into contact with the detector surface is 
determined. The associated measured value is then calculated by interpolating the measured 
values of adjacent beams. The measiared value which can be allocated to the beam of the 
projection that passes the voxel, or the corresponding measured value obtained by 
interpolation, is accumulated on the voxel V(x) in step 205. In step 207 a check is made as to 
whether all projections have been considered. If this is not the case, then the flowchart 
branches to step 203. Otherwise a check is made in step 209 as to whether all voxels V(x) in 
the FOV have been passed through. If this is not the case, then the method continues with 
step 201 . If, on the other hand, all voxels F(x) in the FOV have been passed through, then 
the absorption in the entire FOV is determined and this reconstruction method is terminated. 

In order to determine within the examination area the spatial area in which the 
heart is located, in the following step 1 13 the heart is segmented in the reconstructed three- 
dimensional data record. 

A simple possibility for the segmentation is manual segmentation. A user, for 
example a physician, places markings on the surface of the heart in the three-dimensional 
data record. These markings are connected by lines and thus form a network that represents 
the surface of the heart. 

A fiirther possibility for the segmentation consists in the use of a deformable 
model of the heart which is moved, rotated and scaled in the data record such that the 
correlation between the data record and the model of the heart is maximized. This known 
segmentation method is explained for example in "Deformable models in medical image 
analysis: A survey". Medical Image Analysis, 1(2): 91-108, 1996, to which reference is 
hereby made. 

Another known segmentation possibility is based on a region growing or 
region expansion process, in which a user predefines a so-called seed voxel in an object that 



wo 2004/089217 PCT/IB2004/001061 

10 

is to be segmented, that is to say in this case in the heart. Neighboring voxels of the seed 
voxel are then examined using an association criterion to determine whether they do or do 
not belong to the heart. This association criterion may be for example the fact of whether or 
not they are included within a range of values of the data values of the voxel. If a data value 
5 lies within the range of values then the corresponding voxel is assigned to the heart. In the 
next step, the neighboring voxels of the voxel which has been newly assigned to the heart are 
examined with respect to the association criterion and in some circumstances also Eissigned to 
the heart. This method is repeated until no more neighboring voxels can be assigned to the 
heart. 

1 0 The segmented heart shows the spatial area taken up by the heart in the 

examination area. 

This determination of the spatial area taken up by the periodically moving 
object in the examination area represents only one embodiment. Any method which makes it 
possible to determine this spatial area may be used according to the invention. For instance, 

15 the examination area could also be reconstructed using other known reconstruction 

techniques, such as the n-PI method which is described in "The n-PI Method for Helical 
Cone-Beam CT", R. Proska, Th. K5hler, M. Grass, J. Timmer, IEEE Transactions on 
Medical Imaging, Vol. 19, 848-863, September 2000. The reconstruction coxild also be 
carried out using the method from "A General Cone-Beam Reconstruction Algorithm", G. 

20 Wang, T.H. Lin, P.C. Cheng, D.M. Shinozaki, IEEE Transactions on Medical Imaging, Vol. 
12, 486-496, March 1993. The use of the reconstruction method known as ASSR (Advanced 
Single-Slice Rebinning) would also be possible, which although leading to poorer-quality 
images compared to the abovementioned reconstruction does require less computational 
complexity. This method is published for example in "Advanced Single-Slice Rebinning in 

25 Cone-Beam Spiral CT", M. KachelrieU, S. Schaller, W.A. Kalender, Medical Physics, Vol. 
27, No. 4, 754-772, 2000. A reconstruction method for a circular trajectory, along which the 
radiation source moves relative to the object, is shown for example in "3D Cone-Beam CT 
Reconstmction for Circular Trajectories", M, Grass, Th. Kohler, R. Proska, Physics in 
Medicine and Biology, Vol. 45, No. 2, 329-347, 2000. In reconstmction to determine the 

30 spatial area taken up by the moving object, it is important that the resolution of the 

reconstmcted three-dimensional data record is selected such that the spatial contents can be 
determined for example by means of known segmentation methods. This condition is met, for 
example, when a volume of 20 x 20 x 20 cm^ is represented by 64^ or more voxels. Any 
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method capable of segmenting an object in a three-dimensional data record can be used 
according to the invention for the segmentation. 

Once the heart has been segmented in the three-dimensional data record, the 
geometric center of the heart can be determined in step 1 15 by simple geometric 
5 considerations. Only the geometric center of a known three-dimensional object is to be 

detencnined. The geometric center would be, for example, the center of gravity of this object 
if a spatially constant density of the object were to be assumed. 

In step 1 17, a point in time is allocated to each projection determined in step 
105. For this, it is determined which beam fan of a projection irradiates the geometric center 
10 of the heart. The point in time at which the measured values corresponding to this beam fan 
were acquired is allocated to the respective projection. 

If in other embodiments no parallel rebinning according to step 105 has been 
carried out during the determination of the spatial area taken up by the periodically moving 
object in the examination area, then said parallel rebinning must be carried out prior to step 
15 117. 

For the reconstmction of the periodically moving object while taking into 
account the movement recorded during the measurement, repeating time ranges are selected 
in step 119, where the object has assumed at least a substantially identical object state in each 
time range. The at least one substantially identical object state thus occurs in each time range. 

20 Time ranges in which the object states assumed by the object during these time ranges differ 
from one another as little as possible are preferably selected. For example, two object states 
differ little from one another when the difference between two measured signals that 
characterize the respective object states, that is to say for example in the case of the heart the 
difference in the signals from the electrocardiograph, is small compared to the difference 

25 between the maximum and minimum value of the measured signal detected during a 

measurement. A difference between two measured signals is small, for example, if it is less 
than 1%, 2% or 5% of the difference between the maximum and minimum value of the 
measured signal detected during a measurement. In the reconstruction explained further 
below, use is only made of measured values which were acquired during these time ranges. 

30 In the case of the heart it is usefiil to select a time range of the diastolic phase of the 

movement of the heart, since in said phase the movement of the heart is considerably smaller 
than in the systolic phase. This is shown in figure 5, The period of the electrocardiogram H 
consists of a range Hi with relatively little movement and of a range H2 with a lot of 
movement. In this example of embodiment the time range Hi is selected for the ftirther 
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reconstruction since in this range the object states differ less from one another than in the 
range H2. 

In step 121, those projections which were parallel rebinned in step 105 and 
which were acquired during the time ranges Hi are determined. This is again shown in figure 
5 5. In step 117, points in time to and ti were allocated to the projections Po and Pi. The points 
in time to and ti lie in the time range Hi, so that the projections Po and Pi are used for the 
subsequent reconstruction. By contrast, projections to which the points in time ta and ta were 
allocated are not taken into accoxmt since the points in time tz and ts lie in the time range H2. 
If in other embodiments the measured values have not been multiplied 

10 according to step 107 by the cosine of the cone angle of the beam corresponding to the 

respective measured value during the determination of the spatial area taken up by the object 
in the examination area, then this can be carried out prior to the subsequent back-projection. 
The same applies to the filtering according to step 109. 

The measured values of the projections determined in step 121 are then used in 

15 step 123 for the reconstruction of the absorption distribution in the heart by means of a back- 
projection. The individual steps of this reconstruction are shown in figure 6. 

For this, in step 301 a voxel V(x) is determined within a predefinable field of 
view (FOV). The amoxmt of voxels may be for example 512^ and the FOV may be 20 x 20 x 
20 cm^. In step 303, a projection, that is to say a group of beam fans, is then selected from the 

20 projections determined in step 121, which projection has not yet been used to reconstruct the 
voxel V(x) . If no beam of the projection runs centrally through the voxel V(x) , then the 
point at which a central beam would have come into contact with the detector surface is 
determined. The associated measured value is then calculated by interpolating the measured 
values of adjacent beams. The measured value which can be allocated to the beam of the 

25 projection that passes the voxel, or the corresponding measured value obtained by 

interpolation, is multiplied in step 305 by a weighting factor Wp(x) . The subscript index in 
this case denotes the projection selected in step 303. 

The following consideration is used to determine the weighting factor (x) . 
During the acquisition, the voxel V(x) is irradiated for the first time by the cone-shaped 

30 beam cluster at a specific point in time. This point in time is referred to as the sunrise SR. 
The beam leaves the beam cluster again at a second point in time. This point in time is 
referred to as the sunset SS. A factor (x) is then allocated to the projection selected in 
step 303. This factor is smaller for projections whose point in time allocated in step 117 lies 
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relatively close to SR or SS than for projections whose corresponding point in time does not 
lie relatively close to SR or SS. Relatively close may mean for example that the point in time 
of the projection is located in a range that adjoins SR or SS, which corresponds to 5%, 10%, 
15% or 20% of the overall range between SR and SS. One exemplary profile of the factor 
(x) as a function of the position of the point in time of the projection between srmrise and 

sxmset is shown in figure 7. In said figure, the time range between SR and SS is divided into 
three sections. In the first section Ti, which takes up 10% of the time range between SR and 
SS, the factor a^(x) increases from SR up to a value . In the following range T2, the 

factor is constantly equal to and in the third range T3, which also takes up 10% of the 

time range between SR and SS, the factor decreases again toward SS to zero. 

It is furthermore assumed that the movement of the heart is less pronoimced at 
a point in time which lies in the center of the time range Hi, which comprises little 
movement, than at a point in time lying closer to the boundary with the time range H2, which 
comprises a lot of movement Therefore, measured values of projections whose points in time 
lie in the center of the range Hi are more highly weighted than measured values of 
projections whose points in time lie closer to the boundary with the time range H2. In order to 
take this into account, a second factor (x) is assigned to the measured value. The 

subscript index p in this case denotes a projection and, since a point in time is allocated to 
each projection, also a point in time. Depending on the position of the point in time within 
the time range Hi, this factor may have a profile that ensures the following: the factor is 
greater for measured values of projections whose points in time lie in the center of the time 
range Hi than for measiored values of projections whose points in time lie closer to the 
boundary with the time range H2. An exemplary profile is shown in figtire 8. 

Since the same point in time has been allocated in step 117 to all measured 
values whose beams belong to a projection, all measured values of a projection have the same 
weighting factor (x) for a voxel K(x) , which weighting factor is defined by the following 

equation: 

Here, ^ a, (x) JS^ (x) denotes a sum of all projections which have been determined in step 
121 and are not redundant. 
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In step 307, the weighted measured value is accumulated on the voxel V(x) . 
In step 309 a check is made as to whether all projections have been considered. If this is not 
the case, then the flowchart branches to step 303. Otherwise a check is made m step 3 11 as to 
whether all voxels F(x) in the FOV have been passed through. If this is not the case, then the 

5 method continues with step 301 . If, on the other hand, all voxels F(x) in the FOV have been 
passed through, then the absorption in the entire FOV is determined and the reconstruction 
method is terminated (step 313). 

In other embodiments, step 123 and steps 301 to 3 13 may be replaced by other 
known reconstruction methods which generate a three-dimensional data record from the 

10 parallel projections determined in step 121 . 

Furthermore, the reconstruction in steps 301 to 3 13 may be restricted to 
projections which are not redundant. A projection is redimdant with respect to a voxel K(x) 
if the beam of the projection which runs through this voxel V(x) is oriented m the opposite 
direction with respect to a beam of another projection that has already been used to 

1 5 reconstmct this voxel. 



